The First High-Quality Genome Assembly of Freshwater Pearl Mussel Sinohyriopsis cumingii: New Insights into Pearl Biomineralization.

China leads the world in freshwater pearl production, an industry in which the triangle sail mussel (Sinohyriopsis cumingii) plays a pivotal role. In this paper, we report a high-quality chromosome-level genome assembly of S. cumingii with a size of 2.90 Gb-the largest yet reported among bivalves-and 89.92% anchorage onto 19 linkage groups. The assembled genome has 37,696 protein-coding genes and 50.86% repeat elements. A comparative genomic analysis revealed expansions of 752 gene families, mostly associated with biomineralization, and 237 genes under strong positive selection. Notably, the fibrillin gene family exhibited gene family expansion and positive selection simultaneously, and it also exhibited multiple high expressions after mantle implantation by transcriptome analysis. Furthermore, RNA silencing and an in vitro calcium carbonate crystallization assay highlighted the pivotal role played by one fibrillin gene in calcium carbonate deposition and aragonite transformation. This study provides a valuable genomic resource and offers new insights into the mechanism of pearl biomineralization.

A Novel Kunitz-Type Serine Protease Inhibitor (HcKuSPI) is Involved in Antibacterial Defense in Innate Immunity and Participates in Shell Formation of Hyriopsis cumingii.

Serine protease inhibitors (SPIs) are abundantly reported for its inhibition against specific proteases involved in the immune responses, but SPI data related to calcareous shells are scarce. Previously, our research group has reported the proteome analysis of non-nucleated pearl powder, and a candidate matrix protein containing two Kunitz domains in the acid soluble fraction caught our attention. In the present study, the full-length cDNA sequence of HcKuSPI was obtained from Hyriopsis cumingii. HcKuSPI was specifically expressed in the mantle, with hybridization signals mainly concentrated to dorsal epithelial cells at the mantle edge and weak signals at the mantle pallium, suggesting HcKuSPI was involved in shell formation. HcKuSPI expression in the mantle was upregulated after Aeromonas hydrophila and Staphylococcus aureus challenge to extrapallial fluids (EPFs). A glutathione S transferase (GST)-HcKuSPI recombinant protein showed strong inhibitory activity against the proteases, trypsin and chymotrypsin. Moreover, HcKuSPI expression in an experimental group was significantly higher when compared with a control group during pellicle growth and crystal deposition in shell regeneration processes, while the organic shell framework of newborn prisms and nacre tablets was completely destroyed after HcKuSPI RNA interference (RNAi). Therefore, HcKuSPI secreted by the mantle may effectively neutralize excess proteases and bacterial proteases in the EPF during bacterial infection and could prevent matrix protein extracellular degradation by suppressing protease proteolytic activity, thereby ensuring a smooth shell biomineralization. In addition, GST-HcKuSPI was also crucial for crystal morphology regulation. These results have important implications for our understanding of the potential roles of SPIs during shell biomineralization.

Analysis of the immune function of Caspase-3 in Cristaria plicata.

Caspase-3 is generally considered to be the most important terminal shear enzyme in the process of apoptosis, as well as an important part of cytotoxic T lymphocytes (CTL) killing mechanism, which is confirmed to play an important role in vertebrate cell apoptosis and immune system, and is poorly reported in invertebrates. In this paper, we used bioinformatics to perform amino acid multiple sequence alignment and protein structural domain analysis, and constructed a phylogenetic tree to identify the full-length cDNA of the cloned caspase-3 of Cristaria plicata (Named CpCaspase-3). The expression of caspase-1, caspase-7, caspase-8, and caspase-9 was found to be down-regulated by double-stranded RNA interference of CpCaspase-3 in C. plicata. Some degree of disruption of the caspase signaling pathway occurs. The expression of CpCaspase-3 was affected after injection of Lipopolysaccharide (LPS), Peptidoglycan (PGN), polyinosinic-polycytidylic acid (poly(I:C)), and Aeromonas hydrophila. These results were suggested that CpCaspase-3 was involved in the immune response of C. plicata. The wound recovery process of C. plicata was simulated and CpCaspase-3 was found to promote wound recovery. An autophagy inhibition and autophagy activation model of mussels was constructed, where apoptosis and autophagy undergo crosstalk, and inhibition of autophagy induces the onset of apoptosis, and similarly autophagy activation inhibits the process of apoptosis instead. In addition, a recombinant CpCaspase-3-pEGFP-C1 plasmid was constructed for subcellular localization experiments and found that CpCaspase-3 was distributed in both the nucleus and the cytoplasm. This paper aims to unveil the immune mechanism of C. plicata and provide a theoretical basis for the healthy culture of shellfish.

A trait dataset for freshwater mussels of the United States of America.

The United States of America has a diverse collection of freshwater mussels comprising 301 species distributed among 59 genera and two families (Margaritiferidae and Unionidae), each having a unique suite of traits. Mussels are among the most imperilled animals and are critical components of their ecosystems, and successful management, conservation and research requires a cohesive and widely accessible data source. Although trait-based analysis for mussels has increased, only a small proportion of traits reflecting mussel diversity in this region has been collated. Decentralized and non-standardized trait information impedes large-scale analysis. Assembling trait data in a synthetic dataset enables comparison across species and lineages and identification of data gaps. We collated data from the primary literature, books, state and federal reports, theses and dissertations, and museum collections into a centralized dataset covering information on taxonomy, morphology, reproductive ecology and life history, fish hosts, habitats, thermal tolerance, geographic distribution, available genetic information, and conservation status. By collating these traits, we aid researchers in assessing variation in mussel traits and modelling ecosystem change.

Mitochondrial sequence data reveal population structure within Pustulosa pustulosa.

Unionid mussels are among the most imperiled group of organisms in North America, and Pustulosa pustulosa is a freshwater species with a relatively wide latitudinal distribution that extends from southern Ontario, Canada, to Texas, USA. Considerable morphological and geographic variation in the genus Pustulosa (formerly Cyclonaias) has led to uncertainty over species boundaries, and recent studies have suggested revisions to species-level classifications by synonymizing C. aurea, C. houstonensis, C. mortoni, and C. refulgens with C. pustulosa (currently P. pustulosa). Owing to its wide range and shallow phylogenetic differentiation, we analyzed individuals of P. pustulosa using mitochondrial DNA sequence data under a population genetics framework. We included 496 individuals, which were comprised of 166 samples collected during this study and 330 additional sequences retrieved from GenBank. Pairwise ΦST measures based on ND1 data suggested there may be up to five major geographic groups present within P. pustulosa. Genetic differentiation between regions within Texas was higher compared to populations from the Mississippi and Great Lakes populations, which may reflect differences in historical connectivity. Mitochondrial sequence data also revealed varying demographic histories for each major group suggesting each geographic region has also experienced differential population dynamics in the past. Future surveys should consider exploring variation within species after phylogeographic delimitation has been performed. In this study, we begin to address this need for freshwater mussels via the P. pustulosa system.

Expression of bone morphogenetic protein 10 and its role in biomineralization in Hyriopsis cumingii.

Shells and pearls are the products of biomineralization of shellfish after ingesting external mineral ions. Bone morphogenetic proteins (BMPs) play a role in a variety of biological function, and the genes that encode them, are considered important shell-forming genes in mollusks and are associated with shell and pearl formation, embryonic development, and other functions, but bone morphogenetic protein 10 (BMP10) is poorly understood in Hyriopsis cumingii. In this study, we cloned Hc-BMP10 and obtained a 2477 bp full-length sequence encoding 460 amino acids with a conserved TGF-ß structural domain. During the embryonic developmental stages, the cleavage stage had the highest expression of Hc-BMP10, followed by juvenile clams; the expression in the mantle gradually decreased with increasing mussel age. A strong signal was detected on epidermal cells on the mantle edge by in situ hybridization. In both the shell notching and inserting operations of the pearl fragment assay, we found that the expression of Hc-BMP10 increased after the above treatments. RNA interference assays showed that the silencing of Hc-BMP10 resulted in a change in the morphology of the prismatic layer and nacreous layer, with the prismatic layer less closely aligned and the disordered aragonite flakes in the nacreous layer. These findings indicate that Hc-BMP10 is involved in the growth and development of H. cumingii, as well as the formation of shells and pearls. Therefore, this study provides some reference for selecting superior species for growth and pearl breeding of H. cumingii at a molecular level and further investigation of the molecular mechanism for biomineralization of Hc-BMP10.

Study on the Role of Mitophagy Receptor PHB2 in Doubly Uniparental Inheritance of Hyriopsis cumingii.

In bivalves, the heterogeneity of mitochondrial DNA and its unique mode of transmission have been the focus of attention, which is called doubly uniparental inheritance (DUI). Prohibitin-2 (phb2) is a mitochondrial inner membrane protein that is a key mitophagy receptor for parental mitochondrial removal. Hyriopsis cumingii is a freshwater bivalve in China, the full-length cDNA of H. cumingii phb2 (named Hcphb2) is 2917 bp and encodes a total of 300 amino acids, a highly conserved sequence. Hcphb2 was highly expressed in the ovary. In the gonadal tissues of 5- to 8-month-old female mussels, the expression level of Hcphb2 continued to significantly increase. After Hcphb2 siRNA interference in 6-month-old female mussels, the expression of M-COII, a marker gene on M-type mitochondria, showed a considerable increase (p < 0.05). In contrast, the expression of autophagosome formation and maturation-related genes, atg4b, atg5, atg12, and atg16l, in the ATG family genes was significantly decreased (p < 0.01). Subcellular localization showed that Hcphb2 appeared in spermatogonia, spermatocyte, spermatid, and sperm, and its location changes synchronize with the behavior of M-type mitochondria location changes in DUI species. And it was found that miR-184 negatively regulated Hcphb2. The above results suggest that the mitochondrial autophagy receptor gene Hcphb2 may be associated with the degradation of M-type mitochondria in the freshwater mussel. This process requires multiple genes to participate, of which Hcphb2 and autophagy genes are only some of those that may play a role.

JNK regulates the Nrf2/NQO1-ARE pathway against Microcystins-Induced oxidative stress in freshwater mussel Cristaria plicata.

In response to stress, cells can utilize several processes, such as the activation of the Nrf2/Keap1 pathway as a critical regulator of oxidative stress to protect against oxidative damage. C-Jun N-terminal kinase (JNK), a member of the mitogen-activated protein kinase (MAPK) family, is involved in regulating the NF-E2-related nuclear factor 2 (Nrf2)/antioxidant response element (ARE) signaling pathway. NAD(P)H quinone redox enzyme-1 (NQO1), a downstream target gene of the Nrf2 pathway, plays a vital role in removing peroxide and providing resistance to oxidative injury. We found that microcystins (MCs) stimulated CpNrf2 to express and increase anti-oxidative enzyme activities in a previous experiment. In our current study, the full-length cDNAs of JNK and NQO1 from Cristaria plicata (designated CpJNK and CpNQO1) were cloned. The relative levels of CpJNK and CpNQO1 were high in hepatopancreas. Upon MCs induction, the relative level of CpNQO1 was increased, whereas that of CpJNK was decreased significantly. In contrast, CpNrf2 knockdown upregulated the expression of CpJNK mRNA and phosphorylation of CpJNK protein (Cpp-JNK), but inhibited CpNQO1 expression. Additionally, we found that JNK inhibitor SP600125 stimulated expression of CpNQO1 and CpNrf2 upon exposure to MCs, and we further confirmed that CpNrf2 protein combined with the ARE element in CpNQO1 gene promoter in vitro, and increased CpNQO1-ARE-luciferase activity in a CpNrf2-dependent manner. These findings indicated C. plicata effectively alleviated MC-induced oxidative injury through JNK participated in regulating the Nrf2/NQO1-ARE pathway.

Identification and Characterization of C-Mos in Pearl Mussel Hyriopsis cumingii and Its Role in Gonadal Development.

C-Mos, a proto-oncogene, regulates oocyte maturation by activating the classical MAPK pathway in cells. To examine the function of C-Mos in Hyriopsis cumingii, C-Mos was identified in this study. The full-length cDNA of C-Mos was 2213 bp, including 144 bp in the 5' UTR, 923 bp in 3' the UTR, and 1146 bp in the open reading frame (ORF) region. During early gonad development, the expression of C-Mos from 4 to 6 months of age in H. cumingii was significantly higher than that in other months, with the highest expression in 6-month-old H. cumingii, suggesting that C-Mos may be involved in early gonadal development in H. cumingii. Clear hybridization signals were found by in situ hybridization in the oocytes, oocyte nucleus and oogonium, and a small number of hybridization signals were found in the follicular wall of the male gonads. In addition, the C-Mos RNA interference (RNAi) assay results showed that the knockdown of C-Mos caused a down-regulation of ERK and P90rsk. In summary, these results indicate that C-Mos has a crucial part to play in gonadal development in H. cumingii.

A Comparison of Non-Destructive Visceral Swab and Tissue Biopsy Sampling Methods for Genotyping-by-Sequencing in the Freshwater Mussel Fusconaia askewi.

Limiting harm to organisms caused by genetic sampling is an important consideration for rare species, and a number of non-destructive sampling techniques have been developed to address this issue in freshwater mussels. Two methods, visceral swabbing and tissue biopsies, have proven to be effective for DNA sampling, though it is unclear as to which method is preferable for genotyping-by-sequencing (GBS). Tissue biopsies may cause undue stress and damage to organisms, while visceral swabbing potentially reduces the chance of such harm. Our study compared the efficacy of these two DNA sampling methods for generating GBS data for the unionid freshwater mussel, the Texas pigtoe (Fusconaia askewi). Our results find both methods generate quality sequence data, though some considerations are in order. Tissue biopsies produced significantly higher DNA concentrations and larger numbers of reads when compared with swabs, though there was no significant association between starting DNA concentration and number of reads generated. Swabbing produced greater sequence depth (more reads per sequence), while tissue biopsies revealed greater coverage across the genome (at lower sequence depth). Patterns of genomic variation as characterized in principal component analyses were similar regardless of the sampling method, suggesting that the less invasive swabbing is a viable option for producing quality GBS data in these organisms.

Molecular identity crisis: environmental DNA metabarcoding meets traditional taxonomy-assessing biodiversity and freshwater mussel populations (Unionidae) in Alabama.

The use of environmental DNA (eDNA) to assess aquatic biodiversity is a growing field with great potential for monitoring and managing threatened species, like freshwater mussel (Unionidae) populations. Freshwater mussels are globally imperiled and serve essential roles in aquatic systems as a food source and as a natural water filter making their management essential for ecosystem health. Unfortunately, mussel populations are often understudied, and challenges exist to accurately and efficiently describe the full suite of species present. Multispecies eDNA approaches may also be more challenging where freshwater mussel populations are most diverse due to ongoing and significant taxonomic restructuring that has been further complicated by molecular phylogenies using mitochondrial genes. For this study, we developed a microfluidic metabarcoding array that targets a wide range of species, from invertebrates to fishes, with an emphasis on detecting unionid mussels known to be present in the Sipsey River, Alabama. We compared mussel species diversity across six sites with well-studied mussel assemblages using eDNA surveys and traditional quadrat surveys in 2016. We examined how factors such as mussel population density, biomass and location in the river substrate impacted our ability to detect certain species; and investigated unexpected eDNA detections through phylogenetic analysis. Our eDNA results for fish and mussel species were broadly consistent with the data from traditional electrofishing and quadrat-based field surveys, although both community eDNA and conventional sampling detected species unique to that method. Our phylogenetic analysis agreed with other studies that treat Pleurobema decisum and P. chattanoogaense as synonymous species; however, they are still listed as unique species in molecular databases which complicates their identity in a metabarcoding assay. We also found that Fusconaia flava and F. cerina are indistinguishable from one another using a portion of the NADH dehydrogenase Subunit 1 (ND1) marker, which may warrant further investigation into whether or not they are synonymous. Our results show that many factors impacted our ability to detect and correctly identify Unionidae mussel species. Here we describe the obstacles we faced, including the murky phylogeny of Unionidae mussels and turbid river conditions, and our development of a potentially impactful freshwater mussel monitoring eDNA assay.

The thioredoxin expression of Cristaria plicata is regulated by Nrf2/ARE pathway under microcystin stimulation.

Thioredoxin plays an important role in inhibiting apoptosis and protecting cells from oxidative stress. This study was aimed to clarify how the expression of Trx from Cristaria plicata is regulated by Nrf2/ARE pathway. The expression of CpTrx mRNA was significantly up-regulated in gill and kidney tissues under microcystin stress. The Nrf2 gene of Cristaria plicata was identified to possess an auto active domain bit. While CpNrf2 was knocked down by specific small RNA, CpTrx mRNA expression was significantly down-regulated. The promoter of CpTrx gene had high transcriptional activity, and this basic transcriptional activity persisted after ARE element mutation. The region of promoter -206 to +217 bp was a core promoter region and had forward regulatory elements. Gel shift Assay exhibited that the CpTrx promoter could bind to the purified proteins CpNrf2 and CpMafK in vitro. The binding phenomenon disappeared after the ARE element mutation in promoter region. Subcellular localization experiments displayed that fluorescence overlap between CpNrf2 and Trx promoter increased under microcystin toxin stress. These results suggested that Trx expression was regulated by Nrf2/ARE pathway under oxidative stress.

Low-complexity domain-containing protein (LCDP), induced accumulation of calcium carbonate individual crystals to irregular polycrystals in the triangle sail mussel Hyriopsis cumingii.

The nacre layer is composed of sheet-like aragonite crystals, with often loosely arranged polycrystal aragonite sheets which may induce poor mechanical properties in shells. In this study, a full-length low-complexity domain-containing protein (LCDP) cDNA from the triangle sail mussel Hyriopsis cumingii was generated and its role in shell formation investigated. The full-length cDNA was 1058 bp; it had an open reading frame (ORF) of 714 bp encoding 237 amino acids and contained a 20-amino acid signal peptide at the N-terminus and two low-complexity domains. H. cumingii LCDP was not homologous with other species. Tissue expression analyses showed that LCDP was specifically expressed in the mantle. In shell repair assays, significantly higher LCDP expression was observed in the shell repair group from days 12-21 (p < 0.01). After LCDP silencing, aragonite flake shapes in pearl layers became irregular with disordered deposition, while calcium carbonate (CaCO3) crystal surfaces in prismatic layers became rougher and organic matrices between crystals appeared skeletonized, indicating the importance of biomineralization. Our in vitro CaCO3 crystallization assays showed that LCDP induced single crystals to polycrystals, probably via loose arrangement between aragonite flakes. These results provide new insights on freshwater mollusk biomineralization and a theoretical basis for improved pearl quality.

Prx5 of Cristaria plicata has antioxidant function and is regulated by Nrf2/ARE signaling pathway.

Cristaria plicata is one of the more important freshwater pearl bivalves in China, which is susceptible to pathogen infection, and greatly impacts the ability of breeding pearls. Nrf2/ARE signaling pathway and its downstream target gene Prx5 have endogenous antioxidant functions to protect cells from oxidative damage. The full-length cDNA of Prx5 was cloned from C. Plicata, which was 1420 bp, encoding a total of 189 amino acids and had two conserved cysteine residues (Cys78 and Cys179). The amino acid sequence of CpPrx5 was highly similar to Prx5 of other species. Real-time fluorescence quantitative PCR showed that CpPrx5 was distributed in various tissues of mussels, and the highest expression was in hepatopancreas. The expression of CpPrx5 up-regulated in hepatopancreas and gills after LPS, PGN and Poly:I:C stimulation. The recombinant plasmid DE3-PGEX-4T-1-CpPrx5 was expressed in Escherichia coli BL21 and showed antioxidant activity. With the increase of CpPrx5 protein concentration, the superhelical form of DNA was protected. The expression of CpPrx5 was up-regulated after interference CpKeap1 and down-regulated after interference CpNrf2. Gel block assay showed that CpNrf2 and CpMafK proteins blocked CpPrx5 promoter. Subcellular localization showed that CpPrx5 was located in 293T nucleus and cytoplasm and CpMafK was located in 293T nucleus. GST-Pull down verified that CpMafK and CpPrx5 could bind in vitro. These results indicated that Prx5 had antioxidant function and could protects DNA from oxidative damage, and participated in transcriptional regulation by combining with the transcription factor MafK. In addition, MafK could combine with Nrf2 to regulate the downstream target gene Prx5.

Identification of Nrf2/Keap1 pathway and its transcriptional regulation of antioxidant genes after exposure to microcystins in freshwater mussel Cristaria plicata.

Microcystins (MC) are one of the most abundant and widely distributed cyanotoxins in aquatic systems. MC inhibits the functions of protein phosphatase 1 and 2A (PP1/2A), which can seriously affect ecosystem integrity. The NF-E2-related nuclear factor 2 (Nrf2)/Kelch-like epichlorohydrin-related protein-1 (Keap1) signaling pathway protects against oxidative damage by activating phase II detoxification/antioxidant enzymes. Our previous study revealed that MC upregulates the expression and enhances the activities of the antioxidant enzymes by stimulating the CpNrf2 signaling pathway. In the current study, to further clarify the regulatory role of Keap1 in response to MC-induced oxidative stress in shellfish, we cloned the full-length cDNA of Keap1a and Keap1b from Cristaria plicata (designated CpKeap1a and CpKeap1b), which are 2952 and 3710 bp peptides, respectively. The amino acid sequence of CpKeap1a and CpKeap1b contained Tram-track and Bric-a-brac (BTB), Intervening region (IVR), and Double glycine repeat (DGR) domain. Additionally, CpKeap1a contained two cysteine residues analogous to Cys-273 and -288 in zebrafish, but CpKeap1b did not. Moreover, CpKeap1a and -1b formed a homodimer and heterodimer, respectively, and also formed a heterodimer with CpNrf2. In the hepatopancreas, the expression levels of CpKeap1a and -1b were the highest, but MC treatment down-regulated the expression of these proteins. Moreover, the transcription of antioxidant enzymes with antioxidant response element (ARE-driven enzymes), including CpMnSOD, CpCu/ZnSOD, CpTRX, CpPrx, CpSe-GPx, and Cpsigma-GST was upregulated by CpNrf2 in the hepatopancreas. Compared with the MC-induced group, CpKeap1a-siRNA1117 injection significantly increased the transcription of mRNAs for ARE-driven enzymes and Nrf2. CpKeap1a-siRNA1117 also enhanced the activities of antioxidation enzymes. These findings demonstrated that Keap1a negatively regulated the expression of Nrf2 protein and MC-induced oxidative stress response in C. plicata. Therefore, we speculated that CpKeap1a promoted CpNrf2 by recognizing and binding MC. These events then protected molluscs from MC-induced oxidative damage.

Phylogeny and biogeography of Indochinese freshwater mussels in the genus Pilsbryoconcha Simpson, 1900 (Bivalvia: Unionidae) with descriptions of four new species.

The body of knowledge regarding the classification and evolution of freshwater mussels in the family Unionidae (Bivalvia) in Indochina has recently increased. However, the taxonomic revision of all extant taxa in the region is still ongoing. In this study, the genus Pilsbryoconcha was revised based on an integrative analysis of shell morphology, biogeography, and molecular data. Multi-locus phylogeny indicated the availability of eight species within the genus. Four previously recognized species are P. exilis (Lea, 1838), P. schomburgki (Martens, 1860) stat. rev., P. linguaeformis (Morelet, 1875), and P. carinifera (Conrad, 1837), while four other species are described herein as P. acuta sp. nov., P. mekongiana sp. nov., P. kittitati sp. nov., and P. hoikaab sp. nov. In addition, the neotype of P. carinifera is also designated to clarify its long taxonomic ambiguity. Divergent time estimation and historical biogeography analysis revealed that Pilsbryoconcha originated in the area now called the Khorat Plateau around the middle of the Eocene (mean age = 43.12 Mya), before its range was expanded across Indochina through a series of complex geomorphological changes of river systems, which also led to diversification of the genus.

Identification and Functional Analysis of MAPKAPK2 in Hyriopsis cumingii.

MAPKAPK2 (MK2) is an important regulator of the p38 mitogen-activated protein kinase (p38 MAPK) pathway, which is involved in a plethora of cellular processes concluding the development of gamete cells in meiosis and resisting pathogenic bacterial infestation. Hyriopsis cumingii is a significant mussel resource in China and a good material for pearl breeding. To explore the role of MK2 in H. cumingii, MK2 was identified and cloned, whose full-length cDNA was 1568 bp, including 87 bp in 5' UTR, 398 bp in 3' UTR, and 1083 bp in the open reading frame (ORF) region, encoding 360 amino acids. The expression of MK2 was the highest in the gills. Meanwhile, there was a significant difference in the gonads. After Aeromonas hydrophila and Lipopolysaccharide (LPS) infestation, the transcript level of the MK2 was upregulated in the gills. It indicated that MK2 might be involved in the innate immune response of H. cumingii after a pathogenic attack. After quantifying H. cumingii of different ages, it was found that the expression of MK2 was highest at 1 year old. In situ hybridization (ISH) results showed that the blue-purple hybridization signal was very significant in the oocytes and egg membranes of the female gonads of H. cumingii. The expression of MK2 increased gradually at the age of 1 to 5 months and showed a downward trend at the age of 5 to 8 months. It was suggested that MK2 might play an important role in the formation of primitive germ cells in H. cumingii. To sum up, MK2 might not only be involved in the immune response against pathogenic bacterial infection but also might play an important role in the development of the gonads in H. cumingii.

Analysis of protein kinase C (HcPKC) gene expression and single-nucleotide polymorphisms related to inner shell color traits in Hyriopsis cumingii.

BACKGROUND: Protein kinase C (PKC) is a multifunctional serine and PKC can phosphorylate serine residues in the cytoplasmic domain of tyrosinase, thereby regulating the activity of tyrosinase. Activated PKC is bound to the melanosome membrane, and unactivated PKC is free in the cytoplasm of melanocytes. In this study, we study the role of PKC gene in the melanin synthesis pathway and its effect on the color of the nacre of H. cumingii. RESULTS: In this study, a HcPKC gene in H. cumingii was cloned and its effects on melanin synthesis and nacre color were studied. HcPKC was expressed in both purple and white mussels, and the level of mRNA expression was higher in the purple mussels than in white mussels. Strong and specific mRNA signals were detected in the dorsal epithelial cells of the mantle pallial layer, indicating that HcPKC may be involved in nacre formation. After SNP association with inner shell color related traits, according to the principle that 0.25 < PIC < 0.5 is medium polymorphism and PIC < 0.25 is low polymorphism, the A + 332G site on the HcPKC gene was a site of moderate polymorphism, and the other four sites were low polymorphism sex sites. There was strong linkage disequilibrium among the five loci. A haplotype was constructed and it was found that the frequency of T1 (AGGAA)in the white population was significantly higher than that in the purple population (P < 0.05). CONCLUSION: The study found that HcPKC of H. cumingii can be used as a candidate gene related to inner shell color, and some of the SNP sites can be used for molecular-assisted breeding in the spinnaker mussel, providing a reference for cultivating high-quality freshwater pearls.

Complete mitochondrial genome of freshwater pearl mussel Lamellidens marginalis (Lamarck, 1819) and its phylogenetic relation within unionidae family.

BACKGROUND: Freshwater mussels play a key role in ecology and are often considered as ecological indicators. Conversely, these molluscs are one of the most threatened groups due to several anthropogenic factors. Knowledge of phylogenetic diversity would assist in formulating effective management and conservation measures. Lamellidens marginalis is one of the most widely used freshwater mussel for pearl production in India. The genomic resources for investigating its evolutionary relationship within the Unionidae family are lacking. METHODS AND RESULTS: In this study, the f-type mitochondrial genome of L. marginalis was sequenced using the Illumina sequencing platform. The length of the mitochondrial genome was 15,732 bp consisting of 23 tRNAs, 2 rRNAs and 13 protein coding genes. The arrangement of genes was UF1 type and gene overlap was observed between trnG and nad1. Comparative analysis with other Unionidae species showed a high divergence rate in nad6 followed by nad2 atp8 and nad5. The phylogenetic tree supported monophyly of the Unioninae subfamily and L. marginalis (Parreysiinae) formed a sister branch to this subfamily. The divergence time of the Parreysiinae from its most recent common ancestor (MRCA) was placed in the Mesozoic era. CONCLUSION: This information will be useful for the understanding the evolutionary pattern of the species of Parreysiinae subfamily.

Resolving species-level diversity of Beringiana and Sinanodonta mussels (Bivalvia: Unionidae) in the Japanese archipelago using genome-wide data.

Accurate species identification is of primary importance in ecology and evolutionary biology. For a long time, the unionid mussels Beringiana and Sinanodonta have puzzled researchers trying to unravel their diversity because of their poorly discernible morphologies. A recent study conducted species delineation of unionid mussels based on mitochondrial DNA variation, opening up a new avenue to grasp species diversity of the mussels. However, mtDNA-based classification may not align with species boundaries because mtDNA is prone to introgression and incomplete lineage sorting that cause discordance between species affiliation and gene phylogeny. In this study, we evaluated the validity of the mtDNA-based classification of unionid mussels Beringiana and Sinanodonta in Japan using mitochondrial sequence data, double digest restriction site-associated DNA library (ddRAD) sequencing, and morphological data. We found significant inconsistencies in the mitochondrial and nuclear DNA phylogenies, casting doubt on the reliability of the mtDNA-based classification in this group. In addition, nuclear DNA phylogeny revealed that there are at least two unionid lineages hidden in the mtDNA phylogeny. Although molecular dating technique indicates that Beringiana and Sinanodonta diverged >35 million years ago, their shell morphologies are often indistinguishable. Specifically, morphological analyses exhibited the parallel appearance of nearly identical ball-like shell forms in the two genera in Lake Biwa, which further complicates species identification and the morphological evolution of unionid mussels. Our study adds to a growing body of literature that accurate species identification of unionid mussels is difficult when using morphological characters alone. Although mtDNA-based classification is a simple and convenient way to classify unionid mussels, considerable caution is warranted for its application in ecological and evolutionary studies.